Consensus for second-order multi-agent systems with position sampled data
Wang Rusheng1, Gao Lixin1, Chen Wenhai1, †, , Dai Dameng2
Institute of Intelligent Systems and Decision, Wenzhou University, Wenzhou 325035, China
College of Physics & Electronic Information Engineering, Wenzhou University, Wenzhou 325035, China

 

† Corresponding author. E-mail: whchen@wzu.edu.cn

Project supported by the Natural Science Foundation of Zhejiang Province, China (Grant No. LY13F030005) and the National Natural Science Foundation of China (Grant No. 61501331).

Abstract
Abstract

In this paper, the consensus problem with position sampled data for second-order multi-agent systems is investigated. The interaction topology among the agents is depicted by a directed graph. The full-order and reduced-order observers with position sampled data are proposed, by which two kinds of sampled data-based consensus protocols are constructed. With the provided sampled protocols, the consensus convergence analysis of a continuous-time multi-agent system is equivalently transformed into that of a discrete-time system. Then, by using matrix theory and a sampled control analysis method, some sufficient and necessary consensus conditions based on the coupling parameters, spectrum of the Laplacian matrix and sampling period are obtained. While the sampling period tends to zero, our established necessary and sufficient conditions are degenerated to the continuous-time protocol case, which are consistent with the existing result for the continuous-time case. Finally, the effectiveness of our established results is illustrated by a simple simulation example.

1. Introduction

In the last few decades, collective behaviors of the multi-agent systems have been attracting a great deal of researchers because of their numerous extensive applications such as sensor networks, formation control, consensus, containment, etc. (See Refs. [1]–[5]). Among the major areas of coordination control of multi-agent systems, consensus is one of the most significant and fundamental problems, which aims to design the distributed protocol such that the final states of all agents can converge to an agreement (Ref. [1]).

As all we know, the design of a consensus protocol and consensus analysis mainly rely on the agents’ dynamics and their interaction topologies. In Ref. [6], the authors investigated the first-order consensus problem by the stochastic matrix approach, which provided the theoretical explanation of the Vicsek model (Ref. [7]). In Refs. [8]–[13], difference second-order multi-agent consensus problems were discussed, which have attracted increasing interest due to their numerous practical applications. Up to now, numerous interesting results have been established for the multi-agent consensus. Some other relevant research topics have also been addressed, such as event-triggered, time-delay, finite-time, input saturation, and communication constraint (See Refs. [4], [13]–[17]).

Most of the present approaches to solve the consensus problems of multi-agent systems are based on the state information. However, in many practical applications, some state variables such as velocity cannot be obtained directly because of the equipment constraints or cost saving. Thus, it may be a good choice to design the observer to estimate unavailable state variables, by which the control law is constructed. Up to now, the observer-based consensus problem becomes a significant topic in the area of multi-agent systems. In Ref. [8], Hong et al. proposed a distributed consensus protocol for first-order following agents based on the velocity observer to estimate the leader’s unavailable velocity, which was extended to the case that the leader agent has general linear dynamics (See Refs. [18] and [19]). To track the accelerated motion leader, a distributed tracking law based on the acceleration observer was proposed for the second-order following-agents in Ref. [11]. For consensus problems with general linear dynamics, different design approaches have been provided to construct the full-order, reduced-order, and function observer-based consensus protocols in Refs. [20]–[23].

In many practical applications, due to the application of digital sensors and controllers, agents might only communicate with their neighbors periodically rather than continuously, that is, only sampled-data at discrete sampling instants are available for the synthesis of control protocols. Sampled-data control has increasingly attracted study, and some interesting and useful results have been established. By using direct discretization, the sampled-data based consensus protocols were proposed in Ref. [24], and the observer-based consensus laws with sampled-data was proposed under the fixed and Markov-switching topology in Ref. [25]. The results of Refs. [24] and [25] are essentially obtained for second-order discrete-time multi-agent system but not for second-order continuous-time multi-agent systems. In Ref. [26], necessary and sufficient conditions were established to solve the consensus problem with second-order dynamics via sampled control. The distributed consensus protocols only with sampled position data were proposed for second-order agents in Ref. [27], and the distributed consensus protocols only with sampled position and velocity data were proposed in Ref. [28]. Some sampled control with time-varying sampling intervals have been investigated in Refs. [29] and [30]. In Ref. [31], the authors studied the second-order consensus with sampled position and velocity data via pinning control. The sampled-data consensus problems with inherent nonlinear dynamics was investigated in Ref. [32].

Motivated by the aforementioned works, this paper investigates an observer-based consensus problem of second-order multi-agent systems with position sampled-data under directed fixed topology. Since the observer and sampled-data are both taken into account in the design of consensus protocols in this paper, our addressed problem becomes more interesting and difficult. The main contribution of this paper is that two kinds of observer-based consensus protocols with sampled data are proposed for the second-order multi-agent systems. Different to Ref. [25], the observers involved in this paper have continuous-time dynamics, the analysis and design are based on the multi-agent systems itself but not on its discretization system. By using matrix theory and the sampled control analysis method, some sufficient and necessary conditions with respect to the coupling parameters, spectrum of the Laplacian matrix and sampling period are established.

The rest of the paper is organized as follows. In Section 2, some preliminaries on algebraic graph theory and problem statement are introduced. Two kinds of distributed observer-based consensus protocols with sampled-data are provided in Section 3. Following that, a simulation example is presented to illustrate the theoretical result in Section 4. Finally, conclusions are drawn in the last section.

Now, some notations used throughout this paper are given as follows. Rm×n, Cm×n are the sets of m × n real matrices and complex matrices, respectively. ‖ · ‖ denotes the Euclidean norm. The transpose of matrix A is represented by AT. In and 0n represent the identity and the zero matrices of dimension n, respectively. 1nRn is the column vector with all elements equal to one. Re(μ) and Im(μ) expressed the real and imaginary part of the complex number μ, respectively. AB denotes the Kronecker product of matrices A and B. A matrix is said to be Schur stable if all its eigenvalues lie inside a unit circle. A polynomial is said to be stable if all its roots have negative real parts.

2. Problem Formulation and Preliminaries

Consider a group of N identical agents, whose dynamics is modeled by second-order dynamics as follows:

where ri(t) ∈ Rn and vi(t) ∈ Rn are the position and velocity of the i-th agent, respectively. ui(t) ∈ Rn is the control input.

The interaction relationships among N agents can be conveniently modeled by a simple weighted digraph. Let 𝒢 = (𝒱, , 𝒜) be a directed weighted graph of order N, where 𝒱 = {ν1, ν2, …, νN} is a group of nodes and 𝒱 × 𝒱 is the set of edges. The adjacency matrix 𝒜 = [aij]N × N of directed graph 𝒢 is defined as aii = 0, if (νj, νi) ∈ , then aij > 0 and aij = 0 otherwise. The neighbor set of node vi is denoted by 𝒩i = {j|(vi, vj) ∈ }. The Laplacian matrix is signed by L = [Lij]N × N with Lii = ∑jiaij and Lij = −aij for ij. For a directed graph 𝒢, the Laplacian matrix L has at least one 0 eigenvalue and all other nonzero eigenvalues have positive real parts. Moreover, matrix L has a simple eigenvalue 0 if and only if the directed graph 𝒢 has a directed spanning tree (See Ref. [1]). For convenience, let μi, i = 1,2,…,N − 1 be all nonzero eigenvalues of Laplacian matrix L. The interaction topology 𝒢 involved in this paper satisfies the following assumption.

Assumption 1 The interaction topology 𝒢 contains a directed spanning tree.

The multi-agent system is said to achieve consensus if and hold for all i, j = 1,2,…,N. We say that the protocol ui(t) can solve the consensus problem, if the closed-loop system achieves consensus.

Due to the wide use of computer-based control, sampled data-based protocols are considered. Let tk = kT, k = 0,1,2,… be the sampling times, where T > 0 is the sampling period. In many applications, the velocity information cannot be available because the velocity measurements are inaccurate or it is not equipped with constrains in space, cost, and weight. In our problem, it is assumed that only the position information ri(tk) can be used for designing the agent’s controller. The main objective of this paper is to design a distributed consensus protocol ui(t) only based on sampled position information to make the closed-loop feedback system achieve consensus.

Before establishing our main conclusions, some Lemmas are introduced, which will be used in the convergence analysis.

Lemma 1[1] The complex coefficient polynomial s2 + as + b = 0, and all roots lie inside a unit circle if and only if (1 + a + b)λ2 + 2(1 − b)λ + ba + 1 = 0 is stable.

Lemma 2[33] Consider a complex coefficient quadratic polynomial as follows:

where θ1, θ0, δ1, and δ0 are real constants. Then, f(s) is stable if and only if θ1 > 0 and .

Lemma 3 Consider a sampled-data based system

Then, system (3) is asymptotically stable if and only if matrix

is Schur-stable.

Proof Integrating both sides of system (3) from tk to t, one has

from which we can get

While system (3) is asymptotically stable, we know that . Thus, we have , which means discrete-time system (5) is asymptotically stable, that is,

is Schur-stable. On the other hand, while

is Schur stable, we have . It is easy to see

is bounded. Thus, we know that .

3. Main results

In this section, the sampled-data consensus problem of the second-order multi-agent system will be investigated under the fixed directed topology. Full-order and reduced-order observer-based consensus protocols are proposed to achieve consensus, respectively.

3.1. Sampled-data consensus protocol with full-order observer

Under the assumption that the interaction topology 𝒢 is directed, consider a full-order observer of agent i with sampled data

where i and i are the estimation of the position ri and velocity vi, respectively. h1, h2 are the positive parameters, which will be determined later.

By adopting the zero-order holder strategy, a distributed consensus protocol is proposed as follows:

where k1,k2 are designed to be positive constants.

Now, a consensus condition is established with respect to the coupling parameters, sampling period and spectrum of the Laplacian matrix.

Theorem 1 Suppose that Assumption 1 is satisfied. Then, the multi-agent system (1) can achieve consensus via control protocol (7) and the position and velocity estimation errors asymptotical converge to zero if and only if k1, k2, h1, h2, and T satisfy the following conditions.

Proof Let eri(t) = j(t) − ri(t) and evi(t) = j(t) − vi(t) be the position and velocity estimation errors, respectively. Then, it follows from Eqs. (1) and (6) that

According to Lemma 3, it is not difficult to see that the position and velocity estimation errors asymptotically converge to zero if and only if the following matrix is Schur-stable

Let . By combining Eqs. (1), (7), and (10), the dynamics of the closed-loop system can be expressed as

where

By defining , the overall system can be rewritten as

While all systems (10), i = 1,2,…,N are stable, it is not difficult to see that the multi-agent system (1) achieves consensus via control protocol (7) if and only if system (13) reaches consensus.

Noticing that the digraph 𝒢 contains a directed spanning tree, thus 0 is a simple eigenvalue of matrix L, and L1N = 0N. Moreover, there exists a nonnegative vector rRN such that rTL = 0N and . Define J to be the Jordan form associated with the Laplacian matrix L, there exist nonsingular matrices TRN×N, YRN×(N − 1) such that

where = diag (J1,J2…,Js), and Jl is the Jordan block corresponding to μl, in which μl are the nonzero eigenvalues of L, with multiplicity Nl,l = 1,2,…,s and . It is easy to see that T−1 has form

Due to T−1T = I, we have rT1N = 1, rTY = 0, S1 = 0, and SY = IN − 1. Then, we can verify directly that

Let ξ(t) = (T−1I4n)φ(t). Then, equation (13) can be represented as the following system

which can be divided into the following two subsystems: one is

and the other one is

where ξ = [ξ0T,ξ1T]T.

Denote eφi as

Obviously, for all eφi = 0 (i = 1,2,…,N) if and only if φi = φj (i,j = 1,2,…,N), that is, system (13) reaches consensus. Let

Then, we have

by which we can get

From the above analyses, we know that multi-agent system (13) can achieve consensus if and only if , that is, system (17) is stable.

According to Lemma 3, it is not difficult to see that system (17) is stable if and only if matrix

is Schur stable. Notice that

is a block upper triangular matrix, whose diagonal block matrices are

Therefore, we know that the control law (7) can solve the consensus problem if and only if

are Schur-stable, that is, all its eigenvalues are within the unit circle. In the sequel, we will prove that

are Schur-stable under conditions (8) and (9).

According to Eq. (12), by a simple calculation, one has

Then, the following matrix form can be obtained

Denote

If all matrices F1 and F2i (i = 1,2,…,N − 1) are Schur-stable, then

is Schur-stable. Thus, the multi-agent system (1) can achieve consensus via control protocol (7) and the position and velocity estimation errors asymptotically converge to zero if and only if F1 and F2i are Schur-stable.

The characteristic equations |sIF2i| = 0 and |sIF1| = 0 of matrices F2i and F1 have the following form respectively

In fact, let s = (λ + 1)/(λ − 1). Then, equation (21) can be transformed to

From Lemma 1, the roots of Eq. (21) are within the unit circle if and only if the roots of Eq. (23) are in the open left-half plane. Thus, F2i is Schur-stable if and only if all the roots in Eq. (23) have negative real parts.

According to Lemma 2, equation (23) is stable if and only if

Using the same way to analyze Eq. (22), we can obtain that the roots of Eq. (22) are within the unit circle if and only if h1/Th2 > 1/2 and T < 2/h1. Therefore, the second-order multi-agent system (1) can achieve consensus via protocol (7) if and only if equations (8) and (9) are satisfied. The proof is completed.  ◼

Remark 1 According to Eq. (10), the position and velocity estimation errors asymptotically converge to zero if and only if F1 is Schur-stable, from which we know that while T < min{2h1/h2, 2/h1} the position and velocity estimation errors asymptotically converge to zero, while h1 = 0 or h2 = 0, F1 is not Schur-stable. To get large T, small h1, h2 should be chosen. Unfortunately, small h1 and h2 will slow down the convergence rate.

While T → 0, the protocol (7) degenerates to the continuous-time observer-based consensus protocol

According to conditions (8) and (9), it is easy to establish the consensus condition for continuous-time observer-based consensus protocol (25) as follows.

Corollary 1 Suppose that Assumption 1 is satisfied. Then, the multi-agent system (1) can achieve consensus via control protocol (25) and the position and velocity estimation errors asymptotically converge to zero if and only if all positive constants satisfy

Proof According to Theorem 1, as T → 0, condition (8) must be satisfied. Condition (9) can be expressed as

As T → 0, condition (27) degenerates to

from which the consensus condition (26) is obtained.  ◼

Condition (26) is the same as the consensus condition for state consensus protocol

which was established in Ref. [10]. By adopting the full-order state observer (25), it need not add an additional condition to reach consensus. Naturally, it is interesting to probe how large T can be taken under the condition (26), which is addressed in the sequel. For convergence, denote

Let be the minimum positive zero of fi(s). Then, we can obtain the following result.

Theorem 2 Suppose that Assumption 1 is satisfied. Then, there exists a positive constant T* such that for any T < T*, the consensus problem can be solved by protocol (7) and the position and velocity estimation errors asymptotically converge to zero if and only if condition (26) holds. Moreover,

Proof The necessity is obvious. Now, we prove the sufficiency. While condition (26) holds, it is not difficult to see that equation (29) is well-defined. Obviously, condition (8) holds while T < T*. On the other hand, noticing that , we have fi(T) > 0, which implies that condition (9) holds. According to Theorem 1, the consensus problem can be solved by protocol (7) with any T < T*. ◼

Considering the special case that the graph 𝒢 is undirected and connected, the following corollary can be easily obtained via Theorems 1 and 2.

Corollary 2 Suppose that interaction topology 𝒢 is undirected and connected. Then the second-order consensus problem can be solved by protocol (7), and the position and velocity estimation errors asymptotically converge to zero if and only if k1, k2, h1, h2, and T satisfy the following condition

Proof While the graph 𝒢 is undirected and connected, we know that the Laplacian matrix L is symmetric and positive-semi-definite. Furthermore, we know that μi, i = 1,2,…,N − 1 are real and positive, that is, Im(μi) = 0 and Re(μi) > 0. Then, the Corollary 2 can be obtained directly from Theorem 2.  ◼

Remark 2 While condition (26) holds, we can obtain that fi(s) > 0 and fi(+∞) < 0 (i = 1,2,…,N − 1), from which fi(s) (i = 1,2,…,N − 1) has at least a positive zero point. Moreover, we know that

Thus, T* defined by Eq. (29) is positive. According to Theorem 2, while the sampling period T is sufficiently small, the consensus problem can be solved by protocol (7) if condition (26) holds. For a given sampling period T, we can design the gain parameters in protocol (7) under undirected and connected interaction topology. Let μmax be the maximum eigenvalue of L. The gain parameters can be chosen by 0 < k2 < 2/max, 0 < k1 < 2k2/T, 0 < h1 < 2/T, and 0 < h2 < h1/T.

3.2. Sampled-data control protocol with reduced-order observer

In some applications, the agent can have access to its position information continuously, but might communicate with the neighbor agents periodically. In this part, by assumption that the continuous-time position information can be used by the agent itself but only sampled-data can be used by the neighbor agents, a reduced-order observer based consensus protocol with sampled-data is proposed.

Consider a reduced-order observer for agent i

where zi(t) is the protocol, i(t) is the restructured velocity of vi(t), g > 0 is coupling gain to be determined.

Similarly, the consensus protocol with sampled-data is proposed as follows:

where k1k2 are designed positive constants.

Then, the consensus condition can be obtained by the following theorem.

Theorem 3 Suppose that Assumption 1 is satisfied. Then, the multi-agent system (1) can achieve consensus via control protocol (32) and the velocity estimation errors asymptotically converge to zero if and only if k1, k2, g, and T satisfy the following conditions

Proof Let ei(t) = i(t) − vi(t). For t ∈ [tk,tk+1), it follows from Eq. (31) that the dynamical equations for ei(t) can be described as

which means that the velocity estimation errors asymptotically converge to zero.

By the definition of ei(t) and the protocol (32), the closed-loop dynamics of system (1) can be expressed as

Then, let ε = [rT,vT,eT]T. From Eqs. (35) and (36), the following closed-loop dynamic can be rewritten as

where

Following the same line of the previous subsection, it is not difficult to see that the multi-agent system (37) achieves consensus if and only if the following system is stable

where is the same as that in Eq. (17).

According to Lemma 3, system (38) is asymptotically stable if and only if matrices

are Schur-stable. Since

Similarly, we can obtain that the above matrices are Schur-stable if and only if

and 1 − gT are Schur-stable. From the proof of Theorem 1, we know that F2i is Schur-stable if and only if the following condition hold

Moreover, due to |1 − gT| < 1, T < 2/g is easily obtained. Thus, second-order multi-agent system (1) can achieve consensus via protocol (32) if and only if equations (33) and (34) are satisfied. The proof is completed.  ◼

Similarly, the following result can be obtained as Theorem 2, and the proof is omitted.

Corollary 3 Suppose that Assumption 1 is satisfied. Then, there exists a positive constant * such that for any T < *, the consensus problem can be solved by protocol (32) and the position and velocity estimation errors asymptotically converge to zero if and only if condition (26) holds. Moreover,

where are the same as in Theorem 2.

Remark 3 Similarly, while T → 0, the protocol (32) degenerates to the continuous-time observer-based consensus protocol

According to conditions (8) and (9), it is easy to obtain the consensus condition for continuous-time observer-based consensus protocol (40) as follows:

3.3. Consensus protocol with time delay

In some applications, there exists the input time delay which will damage system performance. While time delay is considered in protocols (6) and (7), the consensus protocol has form as follows:

In the sequel, we solve the consensus problem by protocol (42). The notations are the same as the proof of Theorem 1. By following a similar analysis process, the multi-agent system (13) can achieve consensus if the following system is stable

from which we can get

where ABC are defined in Eq. (12). Suppose that τ < T. From Eq. (44), one has

It is not too difficult to see that system (45) is Schur-stable if all matrices

(i = 1,2,…,N − 1) are Schur-stable. Denote

By simple calculation, we know that Gi is Schur-stable if and only if all matrices 1 and 2i (i = 1,2,…,N − 1) are Schur-stable. Similar to Lemma 3, we know that if system (45) is Schur-stable, then . Thus, we can obtain the following result directly.

Theorem 4 Suppose that Assumption 1 is satisfied and τ < T. Then, the multi-agent system (1) can achieve consensus via control protocol (42) and the position and velocity estimation errors asymptotically converge to zero if and only if all matrices 1 and 2i (i = 1,2,…,N − 1) are Schur-stable.

Remark 4 Similar to Theorem 1, we can use the characteristic polynomial approach to discuss the Schur stability of 1 and 2i. Since the characteristic polynomial of 2i has degree 4 and complex coefficients, the established condition will be very complicated. Thus, we do not intend to address this issue. For a reduced-observer case, we can establish a similar result as Theorem 4.

4. Numerical example

In this section, a simple simulation example is proposed to illustrate the theoretical result.

Consider a second-order multi-agent system composing of N = 6 agents. For simplicity, let n = 1. Its topology structure is described by a directed graph 𝒢 as shown in Fig. 1, and its Laplacian matrix L is obtained as

Fig. 1. Directed graph 𝒢 for six agents.

Apparently, interaction topology 𝒢 has a directed spanning tree. By simple calculation, the eigenvalues of L are μ1 = 0, μ2 = 4.7178, μ3,4 = 4.7190 ± 2.5214i, μ5,6 = 4.1721 ± 0.2678i.

First, we consider a full-order observer case. The initial positions and velocities are taken as

and

The initial positions and velocities for the observers are taken as

Take k1 = 1.1, k2 = 0.8, h1 = 0.6, h2 = 0.9. According to Theorem 2, we can get T* = 0.3539s by a simple computation. Take T = 0.2s, which satisfies the conditions (8) and (9). The position, velocity, estimated position error, and estimated velocity error are depicted in Fig. 2, which shows that the second-order multi-agent system can reach consensus via the protocol (7) and the estimated errors converge to zero. Apparently, while T = 1.4s, which cannot satisfy the conditions in Theorem 1. The related error trajectories are depicted in Fig. 3, which shows that the multi-agent system cannot reach consensus. Thus, the simulation result is compatible with our obtained theoretical result.

Fig. 2. Convergence results with k1 = 1.1, k2 = 0.8, h1 = 0.6, h2 = 0.9, and T = 0.2 s.
Fig. 3. Convergence results with k1 = 1.1, k2 = 0.8, h1 = 0.6, h2 = 0.9, and T = 1.4 s.

For the reduced-order observer case, we take

We choose the scalar parameters and the sampling period as k1 = 1.4, k2 = 0.9, g = 4. By Corollary 3, we can get * = 0.3141s. Take T = 0.26s < *. Obviously, conditions (33) and (34) hold in this case. The simulation results are depicted in Fig. 4, from which we know that the multi-agent system can achieve consensus and the estimated error of velocity converges to zero.

Fig. 4. Convergence results with k1 = 1.4, k2 = 0.9, g = 4, and T = 0.26s.

Here, we only consider the fixed interaction topology case. Certainly, while the interaction topology arbitrarily switches within a finite set of directed topologies, the multi-agent system may not achieve consensus by our proposed protocols. But, while the interaction topology switches with a large enough dwell time, the simulation results also show that the multi-agent system can achieve consensus via our proposed protocols. Our established consensus condition may be able to be weakened by using the concept of an average dwell time. The average dwell time approach to discuss the consensus problem can be found in Refs. [22] and [23].

5. Conclusions

In this paper, the second-order consensus problem with position sampled-data has been addressed. By assumption that the interaction topology is fixed and directed, the continuous-time full-order and reduced-order observers are adopted to estimate the agents’ unavailable states, by which two kinds of sampled-data based consensus protocols have been constructed. By using matrix theory and sampled control analysis method, some sufficient and necessary consensus conditions about coupling parameters and sampling period have been obtained. As a special case, we can obtain the necessary and sufficient conditions for the continuous-time consensus protocols directly. Our future works will concentrate on the sampled data-based consensus issue with the time-varying sampling interval, general linear dynamics, and nonlinear dynamics.

Reference
1Ren WCao Y2011Distributed Coordination of Multi-agent NetworksLondonSpringer-Verlag
2Cao J FLing Z HYuan Y FCao C 2014 Chin. Phys. 23 070509
3Cavagna AGiardina IGinelli FMora TPiovani DTavarone RAleksandra M 2014 Phys. Rev. 89 042707
4Zhou FWang Z JFan N J 2015 Chin. Phys. 24 020203
5Bai LChen FLan W Y 2015 Chin. Phys. 24 090206
6Jadbabaie ALin JMorse A S 2003 IEEE T. Automat. Control 48 988
7Vicsek TCzirok AJacob E BCohen ISchochet O 1995 Phys. Rev. Lett. 75 1226
8Hong YHu JGao L 2006 Automatica 42 1177
9Hu JHong Y 2007 Physica 374 853
10Yu WChen GCao M 2010 Automatica 46 1089
11Gao LZhu XChen W 2012 Int. J. Control Autom. 10 931
12Xu XChen SGao L 2015 Neurocomputing 256 297
13Cao JWu ZPeng L 2016 Chin. Phys. 25 058902
14Xiao FWang LChen T 2014 IEEE T. Automat. Control 59 756
15Su HChen M. Z Q 2015 IET Control Theory 9 399
16Su HChen M. Z QChen G 2015 Int. J. Robust Nonlin. 25 2375
17Su HJia GChen M Z Q 2015 J. Franklin 352 3504
18Hong YWang X 2009 J. Syst. Sci. Complex. 22 722
19Gao LTang YChen WZhang H2011Kybernetika47773
20Li ZDuan ZChen GHuang L 2010 IEEE T. Circuits-I 57 213
21Zhang HLewis F LDas A 2011 IEEE T. Automat. Control 56 1948
22Gao LXu BLi JZhang H 2015 IET Control Theory 9 784
23Gao LCui YXu XZhao Y 2015 J. Franklin 352 5173
24Cao YRen W 2010 Int. J. Control 83 506
25Chen WLi X 2014 Int. J. Robust Nonlin. 24 567
26Liu HXie GWang L 2010 Int. J. Robust Nonlin. 20 1706
27Yu WZhou LYu XLiu JLu R 2011 Automatica 47 1496
28Yu WZheng WChen GRen WCao J 2013 IEEE T. Ind. Inform. 9 2137
29Gao YWang L 2010 IET Control Theory 4 945
30Wang SXie D 2012 IET Control Theory 6 893
31Zhou BLiao X 2014 Nonlinear Dynam. 78 555
32Wen GDuan ZYu WChen G 2013 Int. J. Robust Nonlin. 23 602
33Park P CHahn V1993Stability Theory